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Abstract: In the chiral limit of QCD, the pion decay constant F can be extracted from 
lattice gauge theory by means of a coupling to isospin chemical potential. Here we compute 
the leading correction due to finite volume in the e-expansion of chiral perturbation theory. 
A comparison is made to recent Monte Carlo data. 
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1. Introduction 

Recently, much effort has gone into numerical lattice gauge theory computations in QCD 
very close to the chiral limit. Simulations have started to use chiral fermions for the dynam- 
ical quarks of masses close to the physical values [Q, ^. It is thus becoming increasingly 
important to explore the analytical tools available for extracting physical observables in 
this region of almost massless quarks. Here we shall focus on one of the most important 
low-energy constants of QCD, the pion decay constant F. Because we consider the chiral 
limit, and because our purpose is to derive expressions useful for a comparison between 
finite- volume lattice gauge theory simulations and analytical predictions, we phrase our 
analysis in terms of the so-called e-regime of QCD. Roughly speaking, we will here perform 
an expansion in 1/L, where L is a typical length scale of the given four- volume V, i.e., 
L = fV4. This is a finite- volume scaling regime of almost-massless QCD. 

To extract the pion decay constant F from lattice data one would like to focus on 
an observable that is particularly sensitive to the value of F. An example of such an 
observable has been given in ref. Q. It is a two-point spectral correlation function of 
the Dirac operator when subjected to an (imaginary) isospin chemical potential fi. The 
advantage of this particular spectral function is that it focuses on a null-effect: when the 
isospin chemical potential vanishes, the two-point correlation function has a peak of zero 
width at coincident points. As soon as the imaginary isospin chemical potential is turned 
on, the width becomes finite. Since it depends only the scaling variable ji'^F'^V , this gives 
a direct way to measure F from lattice data. The precise form of the spectral two-point 
function has been computed to leading order from chiral perturbation theory in the e- 
regime for both the quenched theory and for QCD with two nearly massless quarks. More 
recently, it has been shown how these results can also be derived from a chiral Random 
Two-Matrix Theory ||5|. The Random Matrix Theory approach has the advantage of being 
easily generalizable to any number of light quark flavors Nf, and also to the "partially 
quenched" situation in which the chemical potential is coupled only to valence quarks and 
not to the physical u and d quarks of QCD. In addition, it provides all spectral correlation 
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functions, of arbitrary order. In this way even individual eigenvalue distributions can be 
computed, so that one will be able to extract both the chiral condensate S and the pion 
decay constant F by a fit to the distribution of just one single Dirac operator eigenvalue. 
Any sector of fixed topological charge v can be considered. An alternative using ordinary 
baryon chemical potential (which leads to a complex Dirac operator spectrum) has also been 
considered Q. The equivalence between spectral correlation functions of arbitrarily high 
order computed from either the chiral Random Matrix Theories and the Chiral Lagrangian 
to leading order in the e-regime has recently been shown by Basile and Akemann Q. 

When comparing with lattice data it is crucial to be able to estimate the error due to 
restricting the analysis to leading order in the e-expansion of chiral perturbation theory. 
Here we consider the first correction to the effective field theory at an external isospin 
chemical potential. Actually, the calculation is of more general validity, corresponding to 
any vector source = u^T" on SU{N) or U{N). 

In the next section we briefly review the set-up of the e-expansion in chiral perturbation 
theory. In section 3 we present the result of our calculation and make a comparison with 
preliminary lattice data. Section 4 contains our conclusions. 

2. Chiral Perturbation Theory in the e-Regime of QCD 

We are interested in QCD with two light flavors, but the calculation is easily done for 
an arbitrary number of light flavors Nf. We likewise consider isospin chemical in a more 
general sense: We couple all light quarks to a quark (baryon) charge operator B. For two 
light flavors we have B = i-ia^, where /i is isospin chemical potential and is the third 
Pauli matrix. In the more general case we consider Nj to be even, and B = /icr^ 1. In 
fact, all that is important for the calculation that follows is the condition of vanishing trace, 
TrB = 0. Even if we generalize to a vector source that is not traceless, the correction to 
the zero-mode integral is unchanged, and only the constant part of the action is modified. 
The calculation therefore also applies to, for example, ordinary baryon chemical potential. 
In terms of Dirac operators we are thus dealing with two kinds, 

D^ij^:^ ^ [0(A) + ^/x7o]V'f^ = afV?^ 

D_^l;^^^ = []p{A) - ifijoW^^ = iAL"V-^ (2-1) 
and correspondingly two sets of eigenvalues Here lp{A) is the ordinary Dirac operator 

(n) 

and ipj. denotes the eigenfunctions. The method for determining the pion decay constant 
F from these sets of eigenvalues has been explained in refs. ^. That analysis was 
restricted to the tree-level chiral Lagrangian in the e-regime. In order to use it to compute 
F from a lattice simulation, one has to know the size of the leading correction. That is the 
subject of this investigation. 

We are also interested in a situation where the theory under consideration has Nj light 
physical flavors that do not couple to isospin chemical potential, while A'^ light valence 
quarks do couple to it. Computationally, this situation is of much interest since it means 
that the gauge field configurations that need to be used are ordinary ones, without any 
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reference to isospin chemical potential. Only the smallest eigenvalues corresponding to 



valence quark Dirac operators of the kind need to be computed anew. 

Because we consider the chiral limit and assume that chiral symmetry is broken sponta- 
neously, we phrase the analysis in terms of the effective theory of pseudo-Goldstone bosons, 
the chiral Lagrangian. In the e-regime of QCD [|8|, ^] one performs an expansion around 
the zero momentum modes of the Goldstones, taking into account the non-zero momen- 
tum modes in a perturbative manner. Interestingly, the topological charge of gauge field 



configurations v then plays a highly non-trivial role |10|. In a slightly confusing choice of 
terminology this is known as the e-expansion of chiral perturbation theory. 
Let us consider the [Nj -|- A^)-flavor chiral Lagrangian, 



-Tr 



(Vo?7(x))tVo?7(x) + ^(a,C/(x))ta,[/(a 



(2.2) 



where 



(2.3) 



Here we assume that the mass matrix is diagonal and all its elements are taken in the 
e-regime; MTjV ~ 0(1). The N additional replicated flavors can be used to obtain the 
pertinent expression for the partially quenched theory, after embedding A'^ valence quarks 
into these flavors and taking the replica limit ^ 0. We have singled out the zero- 
component of in order to keep the direct connection to isospin chemical potential. As 
explained in the introduction, this restriction is actually immaterial, and a more general 
traceless vector source will lead to identical results. 

Separating the zero-mode from non-zero modes (denoted by \J and ^{x) respectively). 



\J{x) = Ue-x:p{iV2i{x)/F), 
the partition function in a sector with a fixed topological charge v is written 



(2.4) 



Z 



Nf+N 



dU{det Uy exp 



U(Nf+N) 

X exp 

Ju{Nf+N) 
1 



^Tt[MU + MU^ + ^^Tr[U,B][U\B] 



d^x 



+-Tr[{WBU){(\x)B - 2e(x)Be(x) +Be"(x))] 
+Tido^[B + U^BU,C] 



(2.5) 



Note that the integrals are performed over the U{Nf + N) group manifold. The additional 
kinetic term of the singlet non-zero modes is [0] 



A(0 



a 
2N, 



(9^TVe(x))2 + ^(TV^(x)) 



2A^. 



(2.6) 
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with additional constants a and mo that are needed when Nf = 0. Here Nc denotes 
the number of colors, and 1/Nc can be usefully thought of as an expansion parameter. 
The Lagrangian exhibited in ( |2.5| ) includes all terms of 0{e^) that have non- vanishing 
expectation values in the e-expansion, plus a kinetic term of 0{e^) whose expectation 
value vanishes. The expectation value of its square contributes an order term to F.^ 
The counting rules are as follows: ^ ~ 1/L ~ e, ~ e^, B ~ e^.^ The new scaling variable 
F'^ji^V is thus of order unity in this counting, just like MTiV . 

The tree-level two-point correlation function of the ^ fields can be given compactly for 
any value of Nj, and N [|^. In the e-regime, 

{iij{x)iki{y)) = 6u6jkA{x -y) - 5ij5kiG{x - y), (2.7) 

with A(x) and G{x) defined by 



where the indices i,j--- can be taken both in the valence and sea sectors. 

Quenching artifacts appear in G{x) as double poles. One keeps track of the replica 



limit iV ^ in the external indices of eq. (2.7). Recently, the chiral expansion with 
/i 7^ has been considered in the p-regime by Splittorff and Verbaarschot [14|. In that 
case even the propagator matrix is /i-dependent. Here, to 0{e'^) in the e-expansion these 
/^-dependent terms do not contribute. In effect, the propagator we use to compute the 
one-loop correction below is insensitive to fj,. 



3. The Leading Finite- Volume Correction to F 

With the set-up of the previous section it is now straightforward to compute the leading 
(one-loop) correction to F in the e-regime. To this end, we compute the one-loop con- 
tribution to the partition function itself: We bring down the action S and saturate it to 
first non-trivial order in the fluctuation field ^(x). In this way we obtain the leading finite- 
volume corrections in the e-regime. For the term involving the chiral condensate S this was 
done for the full theory in ref. and for the quenched and partially quenched theories 
in refs. [ p^ [l^ . Here we concentrate on the term Tr([B, C/][B, U^). 

The computation is now simple. Working in the replica formalism we can in one sweep 
compute the correction to F in the full theory, the quenched and the partially quenched 

^This term was omitted in an earlier version of this paper. The correct result was first given by Akemann, 
Basile and Lellouch |12[. 

^ As is well known [Il5| , the counting is more involved in the fully quenched theory due to the additional 
terms in £.q{(,). 
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theories^. We find 



F ^ fI^I- ^(A(0) - i I d'xidoA{x)f)^ . (3.1) 

to this order. This calculation shows that 0{e^) effects can be absorbed in the redefinition 
of F in the LO Lagrangian, or the two-matrix theory, analogous to the case of the 1-loop 



correction to E 15]. Another important observation is that the double-pole terms have 
canceled. There is simply no one- loop correction to F in the quenched limit, which can 
be viewed as Nf above. In the partially quenched theory with both valence and 
sea quarks in the e-regime the shift in F is as given in eq. (p.l|); there is no difference 
with the full theory. An immediate consequence of this is that there is to this order no 
dependence on the poorly determined parameters mo and a in the quenched theory. This 



situation is reminiscent of what happens with current correlators in the e-regime ]17|. We 
can intuitively understand the equality between the full theory and the partially quenched 
theory (in the sense defined above) by noting that the pion loop responsible for the one- 
loop correction to this order is insensitive to the value of isospin chemical potential /i. It 
should therefore not matter whether the gauge field configurations have been generated 
with dynamical fermions associated with chemical potential or not. 

As is well known, the integrated massless pion propagator A(0) is ultraviolet divergent. 
In dimensional regularization it actually becomes finite in four dimensions, and it has been 
computed for various finite- volume geometries in ref. |18|. This can be parametrized in 
terms of the so-called "shape coefficients" /3i and /cqo which are functions of the four- 
geometry only, and which are readily computed for any shape following the discussion in 



lH]. The result is 



F^F, = f(i + ^(^ + -W) (3.2) 

for the case V = L^T. The appropriate finite-volume propagator has recently been com- 
puted by Splittorff and Verbaarschot in the p-regime where fj, is of same order as p. 
With the help of that expression one can find the matching between the e and p regimes 
explicitly. 

It is of interest to see the effect of this one-loop shift in F on presently available data 



I IS]. The one calculation which has been performed to date took Nf = 2 and used a 12^ 
volume at a lattice spacing where aF 0.07. Using Pi = 0.1405, the correction factor 
Fl/F = 1.3. This is uncomfortably large. The correction to the condensate is also large. 



Because of the weak dependence on simulation volume, it is difficult to deal with the 
finite- volume corrections by simply pushing to large volume. (The simulations of Ref. B, 



^Here the term "partially quenched" is used in the sense described above: We consider A''^ valence quarks 
coupled to isospin chemical potential, while the Nf physical quarks do not couple to it. 
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the largest-scale Nf = 2 simulations done with fully chiral lattice fermions done to date, are 
only done at volumes which are about twice the size of those of [|l9| .) Precision tests clearly 
require performing simulations at several volumes and observing behavior consistent with 
Eqs. ( |3.2[ ) and Making the simulation volume asymmetric can help: for example, 

Pi is reduced from 0.1405 to 0.0836 for a lattice of size L'^ x 2L, as used by Ref. |Q while 
k()Q is not changed very much (from 0.070 to 0.083). The effects of higher-order terms in 
the chiral expansion will also come in, as 1/V corrections [20|. Neither of the volumes of 
Refs. [||] or are particularly large by today's standards, but because of their extreme 
sensitivity to topology, these kind of e-regime calculations require the use of lattice fermions 
with very good chiral properties. Of course, this problem becomes more acute for larger 
Nf. Here we have considered QCD with only light u and d quarks. 

Finally, we should point out that the present calculation does not directly prove that 
the method proposed in ref. [0] can be applied at one-loop level by just performing the shifts 



( |3.2| ) and (|3.3|). In principle, one would have to perform a full calculation of the needed 
susceptibility and resolvent in the graded or replicated chiral Lagrangian, and then go to the 
cut at imaginary mass [21|. We know of no indication that the result of such a calculation 
to one-loop order in the e-regime should yield a result different from just performing the 
shifts ( |3.2| ) and (^]^), but to our knowledge it has not been shown explicitly. However, 
there is much circumstantial evidence in favor of this being the case. For instance, spectral 
sum rules evaluated to one-loop accuracy are consistent with simply performing such 
a one-loop shift. This holds in the theory with finite chemical potential /x as well [p^]. 



4. Conclusions 



We have computed the one-loop correction to F in the e-regime of QCD with a vector-like 
source that couples like (imaginary) isospin chemical potential. This correction is needed 
in order to quantify the finite-volume correction to the leading-order formula for F in 
the e-regime based on the method described in Ref. Lattice simulations which aim 
to predict low-energy constants from simulations in the e— regime will clearly need to be 
carefully designed to deal with these finite volume corrections. 
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